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The density matrix, i.e. the Fourier transform of the momentum distribution, is obtained an- 
alytically for all magnetization of the Gutzwiller wave function in one dimension with exclusion 
of double occupancy per site. The present result complements the previous analytic derivation of 
the density matrix for the majority spin. The derivation makes use of a determinantal form of the 
squared wave function, and multiple integrals over particle coordinates are performed with the help 
of a diagrammatic representation. In the thermodynamic limit, the density matrix at distance x is 
completely characterized by quantities v c x and v s x, where v s and v c are spin and charge velocities 
in the supersymmetric t-J model for which the Gutzwiller wave function gives the exact ground 
state. The present result then gives the exact density matrix of the t-J model for all densities 
and all magnetization at zero temperature. Discontinuity, slope, and curvature singularities in the 
momentum distribution are identified. The momentum distribution obtained by numerical Fourier 
transform is in excellent agreement with existing result. 

PACS numbers: 71.10.-w, 71.27. +a 



I. INTRODUCTION 



The Gutzwiller wave function is the simplest many-particle wave function capable of capturing both itinerant and 
localized characters of strongly correlated electrons [l| . For a long time, it has been considered as a variational wave 
function for standard models such as the Hubbard model and the t-J model It is now known that the Gutzwiller 

wave function without double occupancy for each site constitutes the exact ground state for the supersymmetric t-J 
model in one dimension, for which both the exchange and transfer terms in the Hamiltonian decay as inverse square 
of the distance [4j. Given this context, it is highly meaningful to derive exact properties for the Gutzwiller wave 
function without double occupation per site. Various quantities have been studied exactly such as the spin and charge 
correlation functions, and the momentum distribution which is the Fourier transform of the density matrix. Although 
the correlation functions have been obtained in a simple closed form I5L l(| , most results for the momentum distribution 
still involve a final integration, or summation over infinite series 0, ID, Q • 

For a case with finite magnetization, Kollar and Vollhardt 0] extended the method of Ref.0], and obtained the 
infinite series expansion of the momentum distribution. The expansion is valid for general dimensions and general 
degree of restricting the double occupation. In one dimension they could sum up the series into a form involving 
a single integral over elliptic functions. The resultant mathematical form is very complicated, and it is difficult to 
obtain any insight from the final formula in Ref.0|. 

On the other hand, a completely different approach has been taken by Arikawa et al 0, |T(|. Namely, the exact 
Green function is integrated over the energy to give the momentum distribution. The result has been obtained only 
for the singlet ground state, and still involves integrals. However, the structure of the integral makes it clear how the 
elementary excitations of spin and charge determine the momentum distribution. 

In a previous paper [111 ], hereafter referred to as I, we adopted another approach, and derived the exact density 
matrix of the Gutzwiller wave function in a closed form. Our results in I, however, have been restricted to the 
majority-spin component. In this paper, we derive the minority spin component, completing the analytic derivation 
for all magnetizations and all densities. In the singlet ground state without magnetization, the limit from the minority 
spin agrees with that from the majority spin. 

Our strategy in I and this paper is to approach in the real space instead of the momentum space. In the real space, 
all singularities in the momentum distribution appear in the asymptotic behavior of the density matrix. Hence we can 
avoid dividing the momentum space into different regions separated by characteristic momenta. This paper gives the 
explicit form of the density matrix with use of Bcsscl functions. The form obtained allows interpretation in terms of 
elementary excitations of spin and charge. It is remarkable that different spin and charge velocities combine in such 
a way as to reproduce the asymptotic behavior expected from combination of the Fermi momenta of up and down 
spins. 



2 



II. GUTZWILLER WAVE FUNCTION IN JASTROW FORM 



We work with the particular representation of the Gutzwiller wave function with exclusion of double occupancy at 
each site. We choose as the reference state the fully polarized state \F) where each site is occupied by an up spin [12j |. 



Then a down spin at site j is created by operating with S*. 
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hj does not obey the ordinary anticommutation rule of fcrmions. We can represent any state in terms of the wave 
function vp({x s },{a; h }) as [13 
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where the set of coordinates {x s } specifies the positions of M magnons, and {x h } specifies those of Q holes. We call 



,t\2 



this scheme the magnon-hole representation. The hard-core constraint is satisfied by the property (Sj) = (hj) 
Sjh)j = 0. By definition ^({a; s }, {x h }) is symmetric against interchange of down-spin coordinates, and antisymmetric 
against hole coordinates. This is related to the commutation rule S~S^ = S~ S~ and the anticommutation rule 



h]h\ 



-h\h\ for i ^ j. Furthermore, we have the relation S i h\ 



hlS~. 
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The Gutzwiller wave function is represented in the magnon-hole representation by 



M+Q 



vMKHz h })= n exp^n^-^n^-^n^-^ 
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where D(i — j) — (L/ir) sin [n (i — j) /L], xj denotes both hole (j — 1, . . . , Q) and magnon (j 



;M + Q) 



coordinates. This form is valid for non-negative magnetization, i.e. 2M + Q < L. We have omitted a normalization 
factor here. Here we have assumed the simplest case of L and M even, and Q odd [l3j], for which the ground state 
is nondegenerate. In the thermodynamic limit, the result does not depend on the choice of even or odd numbers 
of particles. The momentum associated with ^ , g({2 ;S }) l^ 11 }) is t because we have taken M + Q odd. However the 
reference state \F) itself has momentum tt since the occupied set of momentum includes one of the Brillouin zone 
boundary tt, which does not cancel with the rest of occupied momentum. Thus the Gutzwiller state |\&g) nas zero 
momentum because of the cancellation of it in \F) and ^g^x 8 }, {x h }). 

It is convenient to introduce the complex coordinate Zj = exp(2iriXj /L) for both magnons and holes using the 
identity: 



2i sin [tt (xi — Xj) / L] = (zi — Zj)/ ^JZiZj. 



(3) 



Then the Gutzwiller wave function is given as a polynomial of these complex coordinates, apart from the factor for 
the Galilean boost. We work with the form 



vMM) - * m + h ({4)*m({4), 

where the factors in the right-hand side are given by 

Q+M 



(4) 
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i=Q+l Q<i<j<Q+M 
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(0) 



where we have dropped factors of L/ (2tt) from each (zi — Zj ) as compared with Eq. ([2|) . 

The Gutzwiller wave function given by Eq.((2|) turns out to give the exact ground state of the t-J model with a 
special condition for the transfer and the exchange interaction |4[. The t-J model in general is given by 



H 



tJ 



<r=T4 



+ h.c. 
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(7) 
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where Ci a is the annihilation operator of an electron with spin a at site i, is the number operator, and V is the 
projection operator to exclude double occupancy at each site. We assume that the parameters satisfy the following 
condition 

Uj = Jij/2 = tD(i- j)- 2 (8) 

This model is called the supersymmetric t-J model with inverse-square interaction Here the lattice constant is 
taken as the unit of length. Hence the length L of the system gives also the number of lattice sites. 



III. DENSITY MATRIX IN TERMS OF THE WAVE FUNCTION 



We work with the electron density matrix p a (x) for spin up [a —]) and spin down (a =|) electrons, which is the 
Fourier transform of the momentum distribution function n a (k). Thus we are interested in evaluating 

p a (x) = (cl(x + Xj)c a (Xj)) (9) 

in the ground state. Here we have written c a {xj) for Cj a . The left hand side of the equations is independent of Xj 
because of translational invariance. In I, we have derived pt{x) analytically for a magnetization m = nj — n\ > 0. 
The other part Pi(x) requires more elaborate calculation, which is presented in this paper. 

For the down spin part of the density matrix, we have to evaluate the expectation value in the ground state, with 
M + 1 magnons and Q — 1 holes, of 

Pl (x) = (bth x hlb ) = -{b\hlh x b Q ) + 5 X>0 ^±± = -G L (x) + 4,o^^ (10) 

where the hole operators behave as fermionic. We will refer to Gi(x) as the propagator. Ignoring the Kronecker's 5 
in the last expression, the expectation value forces one of the holes in the ket wavefunction to be at x and one of the 
magnons at the zeroth site, while in the bra wavefunction this is reversed. All other holes and magnons have to be 
matched in the bra and ket wavefunctions. 

Using Eq. ((4j) with Eqs.([5]) and ([6|), we use the normalization 

B(M, Q) = J^±... ^Q^ m+h ( Zl . . . ZM+Q )^ m (z Q+1 . . . z M+Q ). (11) 

(Since f m+ i, and VE^ are products of real factors, complex conjugation is not necessary.) Therefore, evaluating the 
expectation value in Eq. (fT0| in the ground state with Q — 1 holes and M + 1 magnons, we have 



Q-lM + l(l- z x ) 2 f dx 2 dxQ-x dx Q+ i dx Q+M 




y m+h (z2 ■ ■ ■ zq-i, z Q+l . . . z M +Q)^ m (z Q+ i . . . z m +q) (12) 



where z x — exp(2nix / L) , and the factor ^m+h^m represents a system with Q — 2 holes and M magnons. Eq. fH 
can be written in terms of the wavefunction for Q holes and M magnons, with the first two holes at the zeroth and 
x sites, as 

1 Q-1M + 1 / dx 2 dxQ-i dxQ + i dx Q+M 



B{M + 1,Q-1) L L J L " L 

Q+M 



(y Z')(Z' l) 

X Yl — ^7= ^m+h( 1 ^ Z x,Z2---ZQ- 1 ,ZQ +1 ...Z M+Q )^l 1 (zQ +1 ...ZM+Q) (13) 



i=Q+l 

which further simplifies to 

U7 B(M+1,Q-1) L L y XQ ' 01K 1 '\J- Q L +1 ZiV^ j V ' 
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In this last expression, the expectation value is evaluated in the ground state with M magnons and Q holes, with hole 
coordinates at Z\ . . . zq and magnon coordinates at zq+i . . . zq + m- The factors of L associated with the Kronecker's 
6 makes them O(l) in the thermodynamic limit. 

In Eq. (fT4")) . there is a factor for each magnon coordinate which can be expanded as 

2 cos <t> — exp[i(j)]/ Zi — Zi exp[— i<f>], (15) 

after defining z x = exp[2i0], or cp = ixxjL. The product of each of these factors can, inside the expectation value of 
Eq. (|14p . be expressed as a sum of diagrams, similar to I. This is pursued further in the next section. 

IV. EVALUATION OF PROPAGATOR IN DIAGRAMS 
A. Basics of the diagram technique 

We first review how the normalization B{M, Q) is evaluated. This discussion is almost identical to that in Ref.[ll[, 
but is reproduced here because it is essential to understanding the diagrammatic expansion. We make extensive use 
of a determinant representation of | M/ g ({^r} ) | 2 following ref.[5(- We introduce a notation 

det(zi ...,z Q ) = det(z p ^)j,=-(q_i)/2...(q-i)/2, (16) 

where the suffix V means the Vandermonde determinant. In the matrix for the determinant, each row has an entry 
(z p ■ ■ ■ Zq) with p ~ — (Q — l)/2in the first row, and p = (Q — l)/2 in the Q-th (last) row. Furthermore, the confluent 
alternant is introduced by 

det(zi . . . z m +q) = det(z? . . . z p Q , z p Q+l ,pz p Q+1 . . . z p M+Q ,pz p M+Q ) p= _ (2M+Q _ 1)/2 _ {2M+Q _ 1)/2 . (17) 

It can be shown that detyi(zi . . . zm+q) has fourth-order zeros {zi — Zj) A if both Zj and Zj are magnon coordinates, 
second-order zeros (z.i — Zj) if one of them is a hole coordinate, and first-order zeros {z; L — Zj) if both and zj are 
hole coordinates. Since in Eq.([2{ and in Eq.(j4]) is real, we obtain up to a real positive factor 

I*g({4)| 2 = *g(M) 2 = (-1) P(M < Q) det(z! . . . z Q ) det( 2l . . . z M+Q ). (18) 

HereP(M,Q) = M(M-l)/2+(M+Q)(M+Q-l)/2 comes from the factors of (-i) inEqs.© and ©. One can verify 
that there is no additional factor of (—1) from the determinants; taking the diagonal term from both determinants, 
which has a coefficient of +1, one obtains the lowest power possible for z\, followed by the next lowest power for z 2 , 
and so on. On the other hand, in the polynomials in Eq.© and (O, since each (zi — Zj) has i < j, this term is obtained 
by taking — Zj from each such factor. Since each such factor in <3>g is repeated in \1/q, we have an overall coefficient 
of +1. The coefficient of this is therefore also +1, and there is no factor of (—1) in going from the polynomials to the 
determinants. Since (-1)" = 1 for any even n, P(M, Q) is equivalent to [Q 2 + Q(2M - l)]/2 = MQ + Q(Q - l)/2. 
Eq. (|lip can now be written as 

5(M,Q) = (-r/W2> J ^l...^M±Qdct(z 1 ,z 2 ...z Q )det(z 1 ...z M+ Q). (19) 

The first determinant is a sum of terms of the form z Pl . . . z V q , and the second determinant is a sum of terms of 
the form zf ... z q Q Q (p Q+1 - qQ+i)z p Q Q ^ +qQ+1 . . . (p Q+M - qQ+M)z Q Q ^M +qQ+M , where {pi . . .p Q } is a permutation of 
-(Q-l)/2...(Q-l)/2, and { qi . . . q Q+M ,p Q+1 . . .p Q+M } is a permutation of -(Q + 2M - l)/2 . . . (Q + 2M - l)/2. 
We adopt the convention that pi > for the magnons. Here p^s and q^s are the momenta of the holes and magnons. 
Integrating the hole coordinates, Z\ . . . zq, we see that pi + = 0, i.e. the p^s and q^s are equal and opposite, both 
covering the range [—(Q — l)/2, (Q — l)/2]. Now integrating the magnon coordinates, zq+i . . . zm+q, we see that 
Pi +qi = for the magnons too. Since the hole coordinates cover [—(Q — l)/2, (Q — l)/2], the magnon p^'s range from 
(Q + l)/2 to (Q + 2M — l)/2, and the minus sign counterpart. 

Diagrammatically, B(M, Q) can be represented as in Fig[lJ where pi + qi = forces all the lines to be horizontal. 
The holes arc represented by the lines connecting black circles (from the Vandermonde determinant) to white circles 
(from the alternant), while the magnons are represented by lines connecting white circles to white circles. Every 
magnon line contributes a factor of Pi — qi, i.e. its horizontal extent, while every hole line contributes a factor of unity. 
Various permutation symmetry factors have to be included to calculate B(M, Q). The final result is 



B(M, Q) = M\Q\{Q + 1)(Q + 3) . . . (Q + 2M - 1). 



(20) 



(a) 



(b) 



FIG. 1: (a) Diagram for evaluating the product of Vandermonde and alternant determinants. The white circles denote the 
exponents in the alternant determinant, integer spaced from — (Q + 2M — 1) /2 to (Q + 2M — 1)/2. The horizontal positions of the 
circles give the exponents, which are equally spaced horizontally. The black circles denote the exponents in the Vandermonde 
determinant, and are slightly separated from the white circles for clarity. A line connecting a black circle to a white one 
represents a hole coordinate, and one connecting a white circle to a white one represents a magnon coordinate. For clarity, 
some of the hole lines have not been shown: every black circle should be connected to the white circle on the opposite side 
at the same level, (b) Schematic of the same diagram, where a vertical projection has been taken and the points are equally 
spaced on the vertical axis. The actual horizontal separation between two points can be calculated by counting how many 
levels down one side and up the other separate them. The dashed line is the magnon hole boundary. 



As a result of this, Eq. tfH)) can be written as 

Q+M 

G l (x) = -A(M ) Q)((L8 XQt0 )(L5 Xl , x ) J] (2 cos <f> - e*+/zi - z^) ) (21) 

i=Q+l 

with 

MM O) = M + 1 Q- 1 B(M,Q) Q(Q-l) (Q + 1)(Q + 3)...(Q + 2M-1) 

{ ,W) L L B(M + 1,Q-1) L 2 Q(Q + 2)...(Q + 2M) ' 1 ' 

Expanding the product of factors in Eq. (l21D . for each magnon coordinate we have a factor of 2cos</i, — exp[i<fi]/ Zi 
or its inverse. Therefore in the diagrams contributing to the expectation value in Eq. (f2"Tj) . the total momentum for 
each magnon is or ±1, while the total momentum for each hole is still zero. Here we have loosely referred these 
integers to "momentum", although physical momentum should be multiplied by 2ir/L. In addition, two of the holes 
— the first and the Q'th ones — have no constraint on their momentum, because of the Kronecker's 6 in coordinate 
space. In Fig. [TJ each magnon line connects two white circles, but can now be horizontal, ascending one step (going 
from left to right) or descending one step. Each hole line connects a black circle to a white circle, and is horizontal. 
However, for two of the holes the black circle can connect to a white circle at any level. We refer to these hereafter as 
"jokers" , because their momenta can assume any value. For clarity, in subsequent figures the black circles are replaced 
by hatched circles for the jokers. Note that if the joker at x = x connects to a line that ascends n levels, it implies 
that the diagram has an overall factor of z™ = exp[2m^>]. On the other hand, for the joker at x = an ascending line 
yields a factor of unity. This is in contrast to the magnon lines, where an ascending line is associated with a factor of 
exp[i(/>], not exp[2i</>]. As discussed at the beginning of this section, each magnon line also contributes a factor equal 
to its total momentum. 

As a warmup exercise, we evaluate ((L8 XQ ^i)(LS XlyX )(2cos4>) M ) in Eq. (f2~lj) . All the magnon lines are horizontal, 
and have the same contribution as in B(M,Q). Therefore either each joker connects to the white circle at its level, 
or the jokers exchange partners. In the first case, the momenta for both the jokers are zero, as in B(M, Q). Instead 
of the sum over xo and xq, we have a factor of L with each Kronecker's S, which has the same effect. In the second 
case, there is an extra factor of — [exp(2i(f))] qQ ~ qi , where the (—1) factor comes from the permutation of terms within 
the determinants. The expectation value is therefore 

2^ £ (23) 
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FIG. 2: Type-I diagrams. The jokers have shaded black circles. The two diagrams shown in the figure are a complementary pair. 
Both jokers are partnered with white circles below the magnon hold boundary. All magnon lines are horizontal or dimerized. 



which sums to 

(2 cos 4>) 

We are now in a position to evaluate the diagrams in Eq . (|2"T|) . We classify the diagrams into four categories, in 
increasing order of complexity, and evaluate each of them. As in Ref. [1 11 ] , we refer to the boundary between magnon 
and hole momenta in B(M, Q) as the magnon- hole boundary, shown by a dashed line in Fig. [IJ even though in a 
general diagram there can be magnon lines that go below this boundary and hole lines that go above this (for the 
holes that are the two jokers). 



1 j sin 2 Q0 1 
Q(Q- 1)1 sin 2 ^ 



B. Type I diagrams 

In these diagrams, the jokers pair with the white circles at their own level, or exchange partners. Therefore neither 
joker connects to a white circle above the magnon hole boundary. All the other holes have horizontal connections. 
This is a generalization of the diagrams just evaluated for ((L5 XliX )(L5 XQi o)) , except that now the magnon lines need 
not be horizontal. Since all the hole lines stay below the magnon hole boundary, the magnon lines have to stay above 
it. As explained in Ref. fTl| , each magnon line is either horizontal or exchanges partners with an adjacent magnon 
line to form a 'dimer'. As a result, the factor of (2cos</>) M changes to F(M,Q), introduced in I. F(M,Q) describes 
the contribution of magnon pairs and dimers, and is defined by the recursion relation 

F(M, Q) = 2 cos<PF(M - 1, Q + 2) - [1 + 1/((Q + 2) 2 - 1)}F(M - 2, Q + 4), (25) 

with the boundary conditions F(Q, Q + 2M) = 1 and F(l, Q + 2M - 2) = 2 cos cf>. From Ea.([2"i]). the contribution to 
the propagator from type I diagrams is then 



Gi(x) = -A(M, Q)F(M, Q) 




(26) 



The Type I diagrams can be divided into two complementary classes: for every diagram where the lines from the 
two jokers are horizontal, there is a a corresponding diagram where their partners are exchanged. This generalizes 
to the subsequent, more complicated, classes of diagrams: diagrams occur in complementary pairs, where the two 
diagrams in a pair differ in that the partners of the jokers are exchanged. The exchange of partners causes a factor 
of — exp[2i(qg — <?i )</>], where qq — q\ depends on the diagram. If <fi — 0, the diagrams cancel each other, and all 
contributions to G\(x) are zero. This is reasonable, since Gt(x) involves destroying a hole and a magnon in the 
ground state, and the same site cannot have both a hole and a magnon. 
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C. Type II diagrams 



In these diagrams, one joker partners with a white circle above the magnon hole boundary. The white circle that 
it would have paired with now has to be partnered. Recalling that, except for the jokers, holes must be linked with 
horizontal lines, this can happen in one of two ways. Either the other joker pairs with the white circle at its level, 
in which case the first joker must be just below the magnon hole boundary, and the magnon line just above the 
boundary moves down one level. Alternatively, we have the same situation with the partners of the jokers exchanged. 
The two possibilities together are the complementary pair of diagrams just discussed, and are shown in Fig. [3] Since 
magnon lines can only have a total momentum of ±1 or 0, if the joker links to a circle I levels above the magnon hole 
boundary, this affects the I magnons immediately above the boundary, as shown in the figure. Beyond the Z'th level, 
the magnon lines behave as before: they are horizontal or dimerized. 

The contribution from these diagrams to the propagator is 

G,fr) = -MM, Q) t ^(M - Q + a) j£t!±! ■ ■ ■ |±|Ef 

Q-l 

exp(2il(j)) exp(— il<p)< 1 — — — - exp ( — 2i(l + p)4>) \ 

P=l 

1 Q-l 

+ exp(il<j>) exp(— 2il(j>)\ 1 — — Yj exp (2i(l +p)(f>) > + <t> <-> — 



p=l 



(27) 



The factors inside the summation are explained as follows. Let us first consider the case when the joker at x — x is 
just below the left edge of the magnon hole boundary, shown in the first panel in Fig. [3l and links to a white circle 
I levels above the magnon hole boundary. The factor of 1/Q is because this joker could have occupied any of the Q 
black circles below the magnon hole boundary, and is now restricted to one position. The factor of F(M — l,Q + 21) 
is because beyond the Z'th level above the magnon hole boundary, the magnon lines are horizontal or dimerizc as in 
Type I diagrams. The contribution from these magnons is as if the magnon hole boundary was moved up by I levels. 
The factor of Q/(Q+ 1) . . . (Q + 21 — 2)/(Q + 2l — 1) is because each magnon line contributes a factor equal to its total 
momentum. The lowest I magnon lines, which are no longer horizontal, have each been reduced in total momentum 
by unity. The factor of exp(2iZ</>) is because the joker at x = x has a total momentum of I, resulting in a factor of 
(z x ) = exp(2il(j>). The factor of exp(—il<j>) comes because I magnon lines each have total momentum —1, from the 
Vandermonde and alternant determinants, so they must each have taken a factor of — exp(— i(j>) from Eq. (|15p . and 
there is an additional (—1)' permutation factor from the determinants because the line from the joker goes over I 
magnon lines it was previously under. 

The next case is the group of diagrams complementary to the ones just considered, shown in the second panel in 
Fig. [3] This group gives rise to the second term inside the first {} brackets in Eq. ([27|) . If the second joker is p levels 
away from the first one, p is summed from 1 to Q — 1, with a weighting factor of 1/(Q — 1) for each. The total 
momentum of the first joker changes by — (I +p) relative to the first diagram from the exchange of partners, resulting 
in a factor of — exp(— 2i(l +p)(j)). 

We also have to consider diagrams similar to these two classes, but in which the first joker is just below the right 
end of the magnon hole boundary. Now the total momentum of the joker is —I, and the I lowest magnon lines are 
ascending rather than descending. For the complementary diagrams, the total momentum of the first joker changes 
by (l+p). 

The remaining diagrams are similar to the four classes considered so far, except that the first joker is the one at 
x = rather than the one at x = x. It is possible to verify for each group separately that this changes cj) to —<\>. All 
these contributions add up to Eq. (|2T|) . It is possible to verify that the cj) <-> —<t> is equivalent to an overall factor of 2. 

The sums over p in Eq. (|27p can be performed, yielding 

G H {x) = -2A(M, Q) ^ F(M -l,Q + 21) {Q - ^ - - } {Q - - _ ^ 

x ( ex P (^)[l - e **f y Sin ^ - + ^ - A (28) 
Q — 1 sin <p 
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FIG. 3: Type-II diagrams. The jokers have shaded black circles. The two diagrams shown in the figure are a complementary 
pair. One joker is partnered with a white circle above the line, I levels above the magnon hole boundary. All magnon lines that 
are not shown are above this, and are horizontal or dimerized. 




FIG. 4: Type-Ill diagrams. Both jokers are partnered with white circles above the line, h and h levels above the magnon hole 
boundary. All magnon lines that are not shown are below or above the [Zi , Z2] block, and are horizontal or dimerized. 



D. Type III diagrams 

Moving on in order of increasing complexity, we consider diagrams in which both jokers link to white circles above 
the magnon hole boundary. The white circles they would have paired with need partners. For Type III diagrams, we 
assume that these partners come from above the magnon hole boundary. This requires that both jokers should be 
immediately below the boundary. (For Type IV diagrams, we will assume that these two white circles pair with each 
other.) By an extension of the argument for Type II diagrams, both the jokers have to be immediately below the 
magnon hole boundary, on either side of it. We assume that the jokers to the left and right pair with white circles l± 
and I2 levels above the magnon hole boundary respectively. As seen in Fig. |4j both white circles immediately above 
the magnon hole boundary connect to circles below the boundary, so that h and I2 have to be no less than 2. There 
is a block of ascending (for I2 > h) or descending (for l\ > 1 2) magnon lines from Z m j n to l maK , where / m i n and / max 
are the lesser and greater of l\ and I2 respectively. All other magnon lines are either horizontal or form dimers. 
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The total contribution to the propagator from Type III diagrams is 

M M 



G m (x) = A(M,Q)J2J2 



1 



Z 1= 2Z 2 =2 



Q(Q-i) 



g 2 (Q + 2Z min ) (Q + 2Z min + 2) . . . (Q + 2/ max - 2) 
Q + 1 (Q + 2/ min - 1)(Q + 2/ min + 1) . . . (Q + 2l max - 1) 

x exp[i(l 2 - h)4>] I exp(2ih4>) |~1 - exp(-2i(Q - 1 + Z x + l 2 )4>)] + exp(-2iZ 2 </>) [l - exp(2i(Q - 1 + h + l 2 )<j>) 



(29) 



In this equation, there is a factor of (—1) from the permutation of the determinants, which cancels the (—1) in front. 
The 1 /Q(Q—1) is because the location of the two jokers is fixed as being immediately below the magnon hole boundary. 
The factors of F come from the dimerized magnon lines below i m i n and above / max . The factors of (Q + const) are 
from the magnon lines: instead of one line of length Q + 1 immediately above the magnon hole boundary, there are 
two lines of length Q that cross the boundary, and instead of Z max — l m in + 1 horizontal magnon lines from the Z m i n to 
the Z max levels, there are Z max — l m - ln ascending or descending lines. The ascending or descending lines also contribute 
l 2 — h factors of — e*^ or l\ — l 2 factors of — exp[— icj>] respectively; the (— l)'s are absorbed in the permutation factor 
at the beginning. The exp(2iZi^>) comes from diagrams shown in Fig. 31 with the x = x joker at the left end of the 
magnon hole boundary. The complementary diagrams have an additional — exp(— 2i(Q — 1 + + 1 2 )4>)- The remaining 
terms are from the same diagrams, but with the x — x joker at the right end of the magnon hole boundary. Eq. (|29p 
simplifies to 



G m (x) = A(M,Q)J2J2 



Q 



M M f t 

F(l min -2,Q+ 2)F(M - l n 

u/ ± 

; 1= 2/ 2 =2 

(Q + 2l min )(Q + 2l miu + 2)...(Q + 2l 



(Q + 2/ min - 1)(Q + 2Un + 1) • ■ . (Q + 2l max - 
x \ 2 cos((7i + l 2 )(j>) - 2 cos((/i + l 2 + 2Q - 2)<j)) 



1) 



(30) 



E. Type IV diagrams 



Finally, we come to diagrams where both jokers connect to white circles above the line, and the white circles they 
would have partnered with (i.e. at their level) pair with each other to form a magnon line. This is shown in Fig. [S] 
The two jokers now no longer have to be immediately below the magnon hole boundary, but can be at any arbitrary 
depth below it. However, in order for their white circles to pair and form a magnon line, both jokers must be at 
the same level, or differ in level by unity. The three possible cases are shown in Fig. [SJ The contributions to the 
propagator from the first case is 



Q/2-1 M M 

G IVa (x) = -A(M,Q) Y, EE( 2cos 

p=0 Z 1 = 12 2 = 1 



x (Q-l-2p) 



Q(Q - 1 



1,Q)F(M- 

'max •■ 

Q + 21 



(Q + 2/ n 



1) . . . (Q + 2U 



1) 



x exp[i(Z 2 - h)4>]< exp(2i(/i +p)4>) 1 - exp(-2i(h +h+Q~ 1)<P) 



+ exp(-2i(/ 2 +p)0) 1 - exp[2i(h + 1 2 + Q - 1)0) 



(31) 



This is similar to Gm(x), so we only discuss the differences, p is the number of levels below the magnon hole boundary 
that the jokers are placed. The factor of 2 cost/) comes from the magnon line below the magnon hole boundary, as 
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(b) (c) 

FIG. 5: Type-IV diagrams. Both jokers are partnered with white circles above the line, h and h levels above the magnon hole 
boundary. The white circles that they leave unpaired connect with each other, forming a magnon line below the magnon hole 
boundary. This line must be horizontal, or ascending or descending by one level. All three cases are shown in the figure. 



does the factor of Q — 1 — 2p. The four terms within the brackets come from diagrams in the first panel of Fig. [5] with 
the x — x joker on the left hand side, the complementary group of diagrams, and the same with the x = x joker on 
the right hand side. 

The contributions from the second and third cases are 



Q/2-2 M M 



G IVb {x) = A(M, Q) £ EE 

p=i 

x (Q-2-2p) 



1 



P =0 ' 

(Q + 2l min ) ...(Q + 21 



1,Q)F(M- 
-2) 



(Q + 2l n 



1) • . . (Q + 2l r! 



1) 



x exp[i(l 2 -h - l)<j>] { exp(2i(l 1 +p+l)<f>) 1 - exp(-2i(h + l 2 + Q - 1)0) 



+ exp(-2i(Z 2 +p)0) l-exp[2i(Zi+Z 2 + Q-l)0) 



(32) 
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and 



Q/2-2 m At 



G IVc (x) = A(M,Q) EE 

p=i 

x (Q-2-2p) 



A_ ^_ 0(0- l/^ 1 
(Q + 2Z min ) . . . (Q + 2/ 



l,Q) J F(M-Z max ,Q + 2Z max ) 
-2) 



(Q + 2/ min - 1) . . . (Q + 2l max - 1) 
x exp[i(l 2 - h + l)4>]lexp(2i(l 1 + p)(f>) 1 - exp(-2i(/i + l 2 + Q - 1)0) 

+ exp(-2i(Z 2 +p+ 1)0) 1 - exp[2i{h + l 2 + Q - 1)0) 



(33) 



closer to the magnon hole boundary from the boundary. The differences 
The sum over p stops at Q/2 — 2, because the lower joker is p + 1 levels 

is an overall minus sign relative to Gjy a which is a permutation factor, 
respectively, because the magnon line below the magnon hole boundary 
horizontal. Finally, inside the brackets, we have exp(2i(Zi + p+ 1)0) in 

p + 1 levels below its partner, i.e. has total momentum l\ -\-p + 1. (The 



Here p is the distance of the joker that is 
between these and Giy a can be explained, 
below the magnon hole boundary. There 
2cos0 is replaced by exp(— i<j>) and exp(i0' 
is descending or ascending instead of being 
Givb because the joker on the left is h 
cxp(— 2i{l 2 +p + l)</>) in Gjy c is similar.) 

It is tedious but straightforward to carry out the sum over p in Eqs. (|3"T1) . ([32l and (|33|) . All three can be added up, 
and yield 

M M ^ 

G IV {x) = —A(M, Q) E E n F( - lniin ~ Q)Fi * M ~ /max ' Q + 2?max) 

i 1= l l 2 =l ^ 

(Q + 2l mia ) . . . (Q + 2Z max - 2) 



(Q + 2/ n 



1) . . . (Q + 2/ n 



x <^ 2 cos((/i + l 2 - 1)0) - 2 cos((Zi + Z 2 + 2Q - 1)0) 



(34) 



V. THERMODYNAMIC LIMIT 



A. Synthesis of all contributions to Gj. 

The expressions we have obtained for the various parts of Gj.(x) simplify in the thermodynamic limit L — > oo. This 
limit makes — > 0, but Q0 and M0 remain finite, and are given by 

Q0 = v c x, (Q + 2M)0 = v s x. (35) 

Here we have fixed Q/L = 1 — n and M/Z = (n — m)/2, and used the velocities w c = 7r(l — n) for charge and 
v s = 7r(l — m) for spin. Elementary excitations with these velocities appear in the supersymmetric t- J model where 
both exchange and transfer decay as the inverse square of the distance [3]. Also, we obtain from Eq. p^|) in the 
thermodynamic limit, 

and 



(Q + 2) . . . (g + 21 - 2) ^ 1 

(Q + 1)(Q + 3) . . . (Q + 2/ - 1) ~* ^(Q + l){Q + 2l-l) 



(37) 
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if Z/L is finite. Finally, from Ref. [Tl(, in the thermodynamic limit 



F(M,Q) = -y/Q(Q + 2M) 



Q + 2M Q 
M g W J o("2"W ~ ( J ^ y ) 



^b(^w s a;)Jo(iw c x) - Jo(^««2c)Yb(~i; c a;) 



(38) 



where J„ and Y„ are Bessel functions of v-th order. Using Eqs. 
namic limit 



(|3"T|) and (f2"rj|) , we obtain in the thermody- 



sm w c a; 



(39) 



Next we proceed to Gjj(x). In Eq. (|2"5|) , replacing the sum over I with an integral with y = (I + Q/2)tfi, and applying 
Eqs. ([3^|) . (f3"8"|) , and (f3"T|) , we obtain in the thermodynamic limit 



G H (x) - -^3? f^' [y (M) j ( y ) e <* - (J <-» Y) 

^ Jv c x/2 1 Z 



g — iv c x/2 giv c x/2 



sm w c x 



u r x 



(40) 



where 3? takes the real part. Here we have used the form for F(M — l,Q + 2l) that is valid if (Q + 21) /L and (M — l)/L 
are both finite in the thermodynamic limit. This is correct over essentially the entire range of I as L — > oo. We see 
that the expressions in Eqs. (|3"9")1 and (|4"0")) are finite in the thermodynamic limit. As given in Eq. (|A2p . we use the 
indefinite integral result 



dye i yZ (y) = [Z (y)+iZ' (y)]ye i y 



where Zq = Jo or Yq is a Bessel function of zeroth order. Then integration of Eq. (|4"0)) results in 



(41) 



2v n 



Gn(x) = — — 



. v s — v c . v s + v c sm v c x 
sm x — sm x 



v c x 



2tt 



Fo(f)Jo(f)-Jo(f)Fo(f) 



smu r icos v r x 



v c x 



v c sm 2 (v c x) 

2~KX 



YoC-f)Ji( V -f)-M V -f)Y 1 ( V -f) 



(42) 



where we have used the relation Jp(y) = —Ji(y) and Y^(y) = — Y\{y) given by Eq. (|A2[) . and the simple Wronskian 
given by Eq. (|Alj) . 

We continue on to evaluation of Gjjj(x) and Gjv(x), which are more complicated. From Eqs. (|3"0)) and (|3^|) . the 
thermodynamic limits of both Gm(x) and Giv(x) are infinite: A(M, Q) is 0(1/ L), the two F functions are both O(L) 
if all their arguments are O(L), as they generically are, and the polynomials in Q in the numerator and denominator 
have an overall 0(l/L 2 ) behavior. The double sum, over li and I2, makes Gm(x) and Giy(x) of O(L). However, 
Cm and Giv almost exactly cancel each other, and their difference is finite in the thermodynamic limit. Comparing 
Eqs. ([50)1 and (f3"4"| . we see that Gm(x) and Giv(x) are equal and opposite, except for four differences: (i) the sums 
over l\ and I2 have different ranges, (ii) the polynomials in Q in the numerator and denominator are different, (iii) 
the arguments of the first F function are different, and (iv) the arguments of the cosines in the brackets arc different. 
Of these, (i) and (ii) cause negligible difference in the thermodynamic limit, but (iii) and (iv) are significant. We 
consider the effects of all four separately. 

First, if l± = 1, a term that exists in Giv(x) but not Gm(x), there is no 0(L) factor from the sum over l\, and 
F(0, Q) = 1 instead of being 0(L). Therefore the contribution of this term to Givix) is 0(1/L) and can be neglected. 
The same is true for I2 — 1. Second, there is a factor of Q/(Q 2 — 1) in Gin instead of l/Q in Giv- This is equal to 
(1/Q)[1 + l/Q 2 + 0(1/ Q 4 )]. The correction from this is Gm(x)/Q 2 , which is 0(1/ L) in the thermodynamic limit 
and can be neglected. 

To account for the difference (iii), we note that the contribution for the sum of Gin and Giv contains the term 



F(M -h,Q + 2h) [F(l 2 -2,(3 + 2)- F(l 2 - 1, Q)} , 



(43) 
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where we have taken Z m ; n = I2 and / max = h without loss of generality. In the thermodynamic limit, the difference 
AF of the two F's is given by 



AF = J — 

4 V v c 



Y (^v c x)J (y 2 ) - (J «- Y) 



v c x 



Y a (±v c x)4(y 2 )-(J^Y) 



(44) 



where y 2 = (Q/2 + l 2 )(f>- Then the whole contribution from (hi) is given by 
, m 1 



G ^ x ) = u? E E [Y Q (\v s x)M yi ) - (J ~ Y)] |r (y 2 ) 
i 1= n 2 =i 

x sRe ifel+a2) (e~ tv " x - e m " x ) 



Jo(^v c x) + v c xJ' Q {^v c x) 



(J~Y) 



(45) 



where yi = (Q + l\/2)(j). Replacing the sums with integrals over y\ and y 2 , we obtain 



G a (x) 



sin(i; c a;) ( 
2x 2 



v s x/2 



dyi 



v c x/2 



dy2 



v c x/2 

Jo{-v c x) + v c xJ' {^v c x) 



Y (^v s x)j ( yi yyi -(./«-> y) 



Y Q (3fy)e i >»-(J~Y) 



(46) 



where 9 takes the imaginary part. Finally, the effect (iv) of the change in the argument of the cosines is to differentiate 
them, since </> — > 0, yielding 



G b (x) 



1 M I 

TT V C X 



AL 2 



- 1 1 

E E [Yohv s x)J (yi) — (J <-> Y)}[Y (y 2 )J (-v c x) - (J <-> Y)]3e^ 1+ ^ (e"*"- + e iv °*) . (47) 



Z 1 = l ! 2 = 1 

As in the case of G a , replacing the sums with integrals over y\ and y 2 , we obtain 



G b (x) 



v c cos(v c x) / ' 



2x 



dyi 



v c x/2 



Hi 



dy2 



v c x/2 



Y (^-x)Jo( yi )e^ — (J <-> Y)] [>o(yx) Jo(y 2 )e^ — (J «-» Y)] . (48) 



Then we obtain the thermodynamic limit of Gih{x) + Giv{x) = G a (x) + Gb{x). 

In Eqs. (j46|) and (|48l) . integration over yi and y 2 can be performed analytically by the use of Eq . (|4"Tj) . The integral 
over 2/2 has the upper limit y\ and the lower limit v c x/2. Surprisingly, the contribution from the lower limit cancels 
out Gu(x) exactly. The final integration over y\ with terms coming from the upper limit y 2 = y± can be carried out 
by using the following formula, which are derived in Appendix: 







Ijj = 


J v c x 


lyy = 


J v c x 


IjY = 


J v c x 






Iyj = 


J v c x 







2 „V 2y l 



dy y-My) [Mv) sin2j/- Jx(y)cos2y] = -3 [( J (y) - i.Ji{y)fy 
' 1 

dyyY Q (y) [Y (y) sm2y~Y 1 (y) cos2y] = -3 [(Y (y) - zY^y)) 2 y 2 e l2y ] 



v B x/2 
v c x/2 
v 3 x/2 
v c x/2 



dy yJo (y) [Y (y) sin 2y - Y 1 (y ) cos 2y] 



[{My) - iJi(y))(Y (y) - iYi{ y )) y 



2 e l2y ] 



dy yY (y) [Mv) sm2y~ J x {y) cos2y] 



Q[(Y (y)-iY 1 (y))(My)-^My))y 2 



v B x/2 
v a x/2 



v s x/2 
v a x/2 



smv s x — sim c i 
2^ 



smi) s i — sim c i 



2tt 



(49) 
(50) 

(51) 
(52) 



Collecting the terms, we find tremendous cancellation. Furthermore, many combinations take the form of the Wron- 
skian Eq. (|Al[) , and can be simplified. The propagator Gi(x) = Gi(x) + Gn(x) + G a (x) + Gb(x) is now given a closed 



14 




FIG. 6: Distance dependence of Gi(x) = —pi(x) + 5 Xj o(n — m)/2 with n = 0.9,0.7,0.5,0.3, and m = 0.1 for all cases. Only 
integer values of x are physically relevant. For comparison, the data obtained by inverse Fourier transform of Kollar-Vollhardt 
results are also shown by points. 



form: 



Giix) 



1 



Anx 2 

1 

-in 
1 

1 

Anx 



[v c x cos(v c x) — sin(-u c ;r)] [v s x cos(v s x) — sm(v s x)} 



Jo(-^v s x)Y (-v c x) - Yq(^v s x)J {^v c x) 



v c sm(v e x)v s sin(i> s :r 
v c sm(v c x) [sbx{v s x) — v s xcos(v s x)] 
v s sin(i> s :r) [sm(v c x) — v c x cos(v c x)} 



Jo{-jV s x)Yi{^v c x) - Y (]^v s x)J 1 {^v c x) 
Jo{^v c x)Yi{^v s x) - Y {^v c x)Ji(^v s x) 



(53) 



This is the main result of the present paper. 

Figure [S] shows Gj(x) for a few cases of n with m fixed as 0.1. In order to test the results, we have also computed 
the inverse Fourier transform of Ref. Q to the real space, as shown by dots. The agreement between the present 
calculation and Ref. [8| is excellent. Although the curves in Fig[6]are shown as continuous, and all pass through the 
origin because Gi(x — ► 0) = from Eq. ((53|) . the coordinate x is only physical for integer x. This point will next be 
discussed in detail. 



B. Comparison between G-\ and Gj 

The up-spin density matrix is related to the propagator G-\{x) = (h^hx) introduced in I as 



/9 T (.T) = (^(x + Xj)c^(Xj)) =S xfi 



1 (n— m) 



G T (x). 



(54) 



From the definition we obtain G-\ (x = 0) = 1 — n, and p-\{x = 0) = (n + m)/2. On the other hand, using Eq. (|10p we 
obtain p^(x = 0) = (n— m)/2 and G^{x = 0) = 0. We note that G a (x = 0) = G a {x — ► 0) for both up and down spin 
despite the fact that only integer values of x are physical. 

In the singlet case m — 0, the up- and down-spin density matrices should coincide. This implies that G|(x;m = 
0) = G T (a;; m = 0) for x ^ 0, but 



G t (x = 0;m = G) ^ G T (x = 0;m = 0), 



(55) 



because of a difference in the Kronecker's 8 terms. Let us see how these relations appear in our result. At v s = n, 
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which happens for m = 0, Gi(x;m = 0) with integer x 0) simplifies to 



«c(-l) X 



4 
sin 



cos(v c x) — 



sin(i> c ;r) 
v c x 



w¥W¥)-*(¥M¥)} 



(56) 



where the superscript in G 



(r) 



signifies that this reduced form is only valid for m = and integer x ^ 0. This indeed 
0. Thus the present calculation and the previous one consistently cover 



agrees with the result in I for Gf with m 
the whole range of magnetization. 

Figure [7] shows the result for G^{x) given by Eq. (f53")) with n = 0.3 and m = 0. The result for non-integer values 
of x is included to clarify the analytic property, especially for small x. The curve for x > 1 tends to the value 
— pj.(0) = — n/2 = —0.15, but finally goes to zero in the limit x — > 0, as discussed for Fig. [6j The period Ax of 
damped oscillation is given by Air ~ 4/n. ~ 13. One might naively expect from the trigonometric function in Eq. (|56j) 
that 



A ? 2tt 4 

Ax = — = - , 

v c 1 — n 



(57) 



which is not the case. This situation becomes clearer in Fig. [8l which shows the result for very dilute density n = 0.1. 
Since the electron correlation is not important except for small x, the result for G±{x) is well fit by the noninteracting 
result: G{x) — — sin(7rnx/2)/(7ra;). However, the hard-core constraint causes deviation from the noninteracting 
behavior near x = 0. The period of damped oscillation is determined by the Fermi momentum kp — 7m/2, which is 
the same as the Fermi velocity in our unit. The Fermi velocity determines the principal oscillation for a general case 
of n such as n = 0.3. 

For comparison between Gj and G| , Fig. [9] shows Gf (x) given by 



G T (x) 



v c cos(7ra;) 



cos(w c a;) 



sm(v c x) 



v c x 



M¥W¥)-M¥)M¥)} 



+ *){v.(f)Mf)-Mf)Mf)} 



(58) 



where cos(7ra;) is used instead of (— l) x used in I, in order to include non-integer values of x. The main difference from 
Gj,(x) is the finite limiting value of 1 — n as x — > 0. There are also oscillations between neighboring integer values 

of x, unlike for G[(x), but these are not physical. We note that, after replacing (— l) x with cos(7rx), G^\x) is the 
same as G-\{x), and therefore has the same non-zero x — > limit and oscillations, emphasizing the limited equivalence 
between G\ r \x) and Gf(ai). 
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G;(X) 




FIG. 8: Distance dependence of Gi(x) in the singlet limit with n = 0.1, m = (solid line and dots). The shape is reasonably 
well fit by — sin(-7rna;/2) / (ttx) (dashed line), which is the case for non-interacting electrons. The inset shows the expanded view 
for small x. 



0.8 r 




FIG. 9: Distance dependence of G-\(x) with n = 0.9,0.7,0.5,0.3, and m = 0.1 for all cases. Only integer values of x are 
physically relevant, and the lines are shown only for easy tracking of Gf (x) at integers. 



C. Asymptotic behavior 



The long-distance behavior of the propagator shows the singularity of the momentum distribution, as discussed in 
I. Using the known asymptotic form of the Bessel functions, and using the Fermi velocities kp f = Tr(n + m)/2 and 
kp[ = 7r(n — m)/2, we have derived the leading behavior of Gj.(x) for large x as 



-1 



- n)(l - m) 



sin k 



fix - 



ni cos kF[X 



(n — m) z 



2tt(1 - n)(l - m)x 32tt 2 (1 - m) 2 (l - n) 2 x 2 



sin kF[X 



sin[(2fci?t + kpi)x] sin[(2fcpf — kpi)x] 



4ir 2 (l-n) 2 x 2 47r 2 (l - m) 2 x 2 

up to 0(l/x 3 ). For reference, our previous result in I for the up spin is quoted [l4[ to 0(l/x 3 ) as 

l + n 2 /[2(l-n) 2 



GT(*)~-Vr=£ 

TTX VI — Tfl 



sin A; 



F-\X 



ni cos kF]X 
2tt(1 — n)(l — m)x 4ir 2 (l — m) 2 x 2 



■ sin k 



f^x - 



sin(feir| -I- 2kFi)x 
47r 2 (l - n) 2 x 2 



(59) 



(60) 



In the case of singlet m = 0, the present result Eq. (f5"9"|) agrees with Eq.(|60p. 
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D. Singularities in momentum distribution 



The momentum distribution n^(k) is obtained from the Fourier transform of Gi(x) as 

L/2 



ni(k) = — (n — m) — G±(x) cos kx, (61) 

x = ~L/2 

An exp(ifcoa;)/a; p+1 term in Gi(x) corresponds to a discontinuity in the p'th derivative of Tij,(fc) at k — fco- Thus 
n|(fc) and its derivatives are discontinuous at k = ±kp^. In addition, the second and higher derivatives of n±(k) are 
discontinuous at \k\ — \kp^ ± 2fepf|. The discontinuity at fc = fc F j^ is derived from Eq. (|59[) as 

A ni {k Fl ) = n x (fc Fi + 0) - n L {k FL - 0) = -^(l - „)(l - m ), (62) 

while the up spin discontinuity is given in I as 

An T (fc FT ) = - (63) 

which agrees with Anj^/cpjJ with m = 0. Note that the majority spin (|) has a larger discontinuity, giving An^kpi) = 
1 in the case of full polarization m = n. This must be the case, since the system becomes the same as free fermions 
without spin. The next singularity is the discontinuity in the slope, which is on top of Ariy{kp^) and is given by 

A ( dn ^A = 2i (64) 

\ dk ) 27r(l-n) 1 /2(l_ m )i/2- 1 > 

Furthermore, the discontinuities in the curvature are obtained from Eq. (|59[) as 



A > cfn L (2k F , + k n ) \ = + (65 ) 



dk 2 J 4tt 2 (1 -n) 3 / 2 

' d»n;(2^-fc f; ) = _ f~n 



dk 2 J ~ 4tt 2 (1 -to) 3 / 2 ' 

In the case of singlet to = 0, Eq. (f6"5)) corresponds to the momentum k = 3kp, while Eq. fB^)) contributes to the 
discontinuity in the curvature at k = kp. On the other hand, we have obtained in I: 

d 2 n T (fc FT +2k Fi )\ _ I 



dk 2 J 47r 2 (l-n) 3 / 2 (l-TO) 1 /2- 

Eqs. (|6T|) and (|63|) are reduced to the same with to — 0. 

Note that there is no singularity in nj(fc) at k = ±(kp-\ — 2kp{), while there is one in rij(fc) at k — ±(kp± — 2kp^ ). 
This curious asymmetry between up and down spins can be understood by appealing to elementary excitations in the 
supersymmetric t- J model 0, 0, EH- Namely, the fermionic excitations are combinations of charge excitations: holons 
and antiholons, and spin excitations: spinons and antispinons. The threshold momentum and the quantum number 
for each gapless excitation is given by [l5|, Gil 

spinon : ± 7rm/2, with spin 1/2, antispinon : ± 7tto, with spin 1, (68) 

holon : ±7m/2, with charge +1, antiholon : ± 7m, with charge —2. (69) 

The simplest is a holon-spinon excitations which makes up a fermionic hole. The excitation gives the characteristic 
momenta, ±fcpj = ±.ir(n + m)/2 and ±fcf J, = ±7r(n — m)/2. 

On the other hand, an electron addition excitation has the charge —1, and must involve an antiholon. The simplest 
combination consists of 

(holon) + (spinon) + (antiholon). (70) 

The threshold momenta are not only ±kp but those which are reduced to ±3fc F in the singlet limit. For example, we 
obtain 



7rn/2 + 7rm/2 + nn = 7m/2 — irm/2 + 2ix(n + m)/2 = kp^ + 2kpp 



(71) 
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FIG. 10: Plot of ny{k) versus k with n = 0.9,0.7,0.5,0.3, and m = 0.1 for all cases. For comparison, the Kollar-Vollhardt 
results [8| are also shown by solid lines. 

Furthermore, the combination (holon)+(spinon)+(antispinon) can also make up a hole excitation with positive 
charge. However the spin can take only the value +1/2, since the antispinon excitation from the magnon condensate 
accompanies the spin flip from down to up, but not the opposite. Because n^{k) — (cj^Cfcf) involves hole excitations 
with spin —1/2 as a result of acting Cfcj to the ground state, the antispinon cannot participate in the process. In 
n|(/c), on the other hand, excitations accompany the spin change +1/2 by c/4. Then a threshold momentum is given 
by 

—Tm/2 — -Km/2 — Trm = Tr(n — m)/2 — 2~n(n + m)/2 = kp± — 2kpp (72) 

Hence there emerges a weak singularity at this momentum and the minus sign counterpart. 

In order to see the global behavior in the momentum space, we numerically carry out the Fourier transform of 
G^(x), and derive the momentum distribution function n^k). The analytic calculation also seems possible, but is 
extremely complicated (l7j . In this paper we are satisfied with numerical comparison with previous results Q. As 
shown in Fig. [101 the agreement between the present and previous results is excellent. It is difficult to identify the 
tiny discontinuities in the curvature at kpi — 2kp] in Fig llOl However, the discontinuity in the slope at fepj. is clearly 
seen and is characterized by Eq. ([M)l . 

VI. CONCLUDING REMARKS 

We have obtained the analytic form of the density matrix of the Gutzwiller wave function with allowance of 
magnetization. Together with the previous results for the majority spin, we have now expressions for both components 
of spins. The density matrix in the thermodynamic limit is characterized by spin and charge velocities, which represent 
elementary excitations in the supersymmetric t- J model with inverse-square exchange and transfer [il. [TH fig. The 
singularities in the momentum distribution are identified completely. The identification through the asymptotic form 
of the density matrix is more effective than previous methods [18j which try to analyze them in the momentum space. 

Although we have dealt with many complicated terms, the final result seems relatively simple, as given by Eq. (|53[) . 
The simplification emerges only in the thermodynamic limit. Similar simplification has also been seen in the exact 
dynamics, where an extremely complicated expression for finite size leads to a final result in terms of simple form 
factors [9|. In the present case, even after the thermodynamic limit is taken, there is considerable further simplification, 
with double integrals involving products of four Bessel functions eventually reducing to Eq. (j55|) . It is tempting to 
speculate that there may be a way to work directly in the thermodynamic limit without bothering about finite size, or 
even to obtain a completely different real space formulation in which the final results are obtained without 'fortuitous' 
cancellations. 
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APPENDIX A: IDENTITIES INVOLVING BESSEL FUNCTIONS 

Here we summarize the properties of Bessel functions, which are used in this paper. Bessel functions have the 
simple Wronskian: 

J (x)Yj(x) - 4(x)Y (x) = — . (Al) 

7TX 

which comes from Abel's identity. Another important property is the differentials: 

Z' (x)=-Z 1 (x), Z[(x) = Z (x)--Z 1 (x), (A2) 

x 

which apply for cases Z v (x) — J„(x), Y L ,(x), and their linear combinations with u = 0,1. With Eq. (|A2[) we obtain 

^ ([Z (x) - tZ 1 (x)]xe^ = A f[ Zo ( x ) + iZ' (x)]xeA = Z (x)e ix . (A3) 

Hence the right-hand side can be integrated. We introduce the notation: u(x) = Jo(x)e' tx ~ U'(x) and v(x) = 
Y (x)e tx = V'(x), where 

U{x) = [J {x) + U' Q {x)] xe ix , V{x) = [Yo(x) + iY^(x)} xe ix . (A4) 
Then we consider the integral: 

K af} = I dx a'(x)/3(x), (A5) 

J a 

where a and (3 denote either U or V. It is obvious that 

K aa = l - [a(bf - a(a) 2 } . (A6) 



Therefore we obtain 



K uu = 1 [U{bf - U{a) 2 } , K vv = l - [V{bf - V{a) 2 } . (A7) 



The other components are manipulated as 



which leads to 



K uv = I dxU'(x)V(x) = U(x)V(x)\ b a - I dxU(x)V'(x) = U(x)V(x)\ b a - K vu , (A8) 



K uv +K vu = U(x)V(x)\ b a . (A9) 



On the other hand, the difference is organized as 



K uv -K vu = f dxJ {x)e tx xe rx [Y (x) + iY{(x)} - (J <-> Y) 

J a 

,6 

= i / dx[J (x)Y^(x) -Y (x)Jo(x)}xe 2lx 

J a 

= I (e 2lb - e 2la ) , (A10) 



20 



where we have used the Wronskian Eg. ljAip . Hence we obtain 



1 



2ib 



e^) + -[U(b)V(b)-U(a)V(a)}, 
e ««) + I [U(b)V(b) - U(a)V(a)\ . 



U I > 



2tt 
-1 




2 ('fc 



r 1/ 



(All) 



By taking a — v c x/2 and b = v s x/2, the imaginary parts of K a p gives the results for I a p quoted in Eqs. (|49|) - (|52l) . 
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